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Great interest has attached to recent D+Au, ^/s — 200 A GeV data at RHIC, obtained with the 
BRAHMS detector. Between pseudorapidities ry = and r; — 3.2 the appropriately defined ratio 
R[DAu/PP], comparing transverse momentum spectra of D+Au to P+P, exhibits a steady decrease 
with rj. This diminuition is examined within a two-stage simulation, the last stage being a purely 
hadronic, reduced energy cascade. The result is an adequate description of the data including the 
so-called Cronin effect. Additionally there is clear evidence for suppression, in the second stage, of 

^^ [ relatively high transverse momentum r; = leading mesons, i. e. the Cronin effect, only near mid 

(^ . rapidity, is appreciably muted by final state interactions. 

o 
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1—5 ; I. INTRODUCTION 

cn , 

04 , It is instructive to begin with a comparison of the minimum-bias experimental distribution dN'^^/drj for charged 

particles at y^ = 200 A GeV and our eventual attempt to describe this. These are shown in Fig(l) together with 
the equivalent PP measurement without any normalising factors imposed. The data shown are from UA5 T| and 
PHOBOS 2]. Clearly charged particle production in D+Au is considerably enhanced relative to PP; but for large ?/ the 
id, D+Au spectrum asymptotically joins PP, both in the data and in our simulation. It is also evident the ij = 3.2 point 
(^ is appreciably suppressed relative to 77 = 0, for D+Au but not for PP. It will become clear that these pseudorapidity 
\^ spectra, which are in fact integrals of the double differential cross-sections {l/2TTp±)[d?N'^'^/dp±dri], are dominated 
^D , by quite small transverse momenta p±. We examine this self-evident thesis in more detail as wc proceed. 
'^ ■ The BRAHMS collaboration '3e\ has focused on the 77 and p^-dcpendence of the ratio 
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where Ncoii is a calculated number of binary NN collisions occurring in minimum-bias D+Au. They consider also the 
combined 77 and p±^ dependence for varying centralities ci^2 of 
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where the denominator is the same ratio for a peripheral setting. The behaviour of the former ratio we contend is 
determined mainly by the low transverse momentum dynamics and by p± distribution in PP collisions, which is for 
us an input to the nucleus-nucleus simulation. There are important dynamic modifications, e. g. the Cronin effect jj], 
but the relation between low and high p± is to a large extent similar to that in PP. The observed behaviour of Rc^p 
with pseudorapidity and centrality is determined by the asymptotic approach of dN/drj in D+Au to that in PP at 
increasing 77, and the resulting diminished variation with 77. 

The code LUCIFER, developed for high energy heavy-ion collisions has previously been applied to both SPS energies 
v/s = (17.2,20) A GeV 5] and to RHIC energies y^ = (56, 130,200) A GeV HQ- We present a brief description 
of the dynamics of this Monte Carlo simulation. Many other simulations of heavy ion collisions exist and these are 
frequently hybrid in nature, using say string models in the initial state J8[ Ifl. p il llll Il3 . Ila ll4L Il5l| together with final 
state hadronic collisions, while some codes are purely partonic [la IrZl llsL Il9l l2(j| in nature. Our approach is closest 
in spirit to that of RQMD and K. Gallmeister, C. Greiner, and Z. Xu and parallel work by W. Gassing 8, 21, 2^. 
Certainly our results seem to parallel those of the latter authors. 

The purpose of describing such high energy collisions without introducing the evidently existing parton nature of 
hadrons, at least for soft processes, was to set a baseline for judging whether deviations from the simulation measured 
in experiments existed and could then signal interesting phenomena. The division between soft and hard processes, the 
latter being in principle described by perturbative QCD, is not necessarily easy to identify in heavy ion data. For the 



relatively simple D+Au system we are interested in separating the effects of our second stage, a lower energy hadronic 
cascade, from those of our first stage, a parallel rather than sequential treatment of initial (target)- (projectile) NN 
interactions. 

II. THE SIMULATION 
A. Stage I 

The first stage I of LUCIFER considers the initial interactions between the separate nucleons in the colliding ions 
A+B, but is not a cascade. The totality of events involving each projectile particle happen essentially together or one 
might say in parallel. Neither energy loss nor creation of transverse momentum {pi_) are permitted in stage I, clearly 
an approximation. A model of NN collisions 5, d|, incorporating most known inclusive cross-section and multiplicity 
data, guides stage I and sets up the initial conditions for stage IL The two body model, clearly an input to our 
simulation, is fitted to the elastic, single diffractive (SD) and non-single diffractive (NSD) aspects J23| of high energy 
PP collisions P, |2J| and PP data [23 . It is precisely the energy dependence of the cross-sections and multiplicities 
of the NN input that led to our successful prediction 0j @] ^^ ^^^ rather small (13%) increase in dN'^'^/drj between 
^/i = 130 and ^/s = 200 A GeV, seen in the PHOBOS data ^. 

A history of the collisions that occur between nucleons as they move along straight lines in stage I is recorded and 
later used to guide determination of multiplicity. Collision driven random walk in p± fixes the p± to be ascribed to 
the baryons at the start of stage II. The overall multiplicity, however, is subject to a modification, based, as we believe 
on natural physical requirements '6|. 

If a sufficiently hard process, for example Drell-Yan production of a lepton pair at large mass occurred, it would 
lead to a prompt energy loss in stage I. Hard quarks and gluons could similarly be entered into the particle lists and 
their parallel progression followed. This has not yet been done. One viewpoint and justification for our approach is to 
say we attempt to ignore the direct effect of colour on the dynamics, projecting out all states of the combined system 
possessing colour. In such a situation there should be a duality between quark-gluon or hadronic treatments. 

The collective/parallel method of treating many NN collisions between the target and projectile is achieved by 
defining a group structure for interacting baryons. This is best illustrated by considering a prototype proton-nucleus 
(P-l-A) collision. A group is defined by spatial contiguity. A proton at some impact parameter b{x±) is imagined to 
collide with a corresponding 'row' of nucleons sufficiently close in the transverse direction to the straight line path 
of the proton, i. e. within a distance corresponding to the NN cross-section. In a nucleus-nucleus (A+B) collision 
this procedure is generalized by making two passes: on the first pass one includes all nucleons from the target which 
come within the given transverse distance of some initial projectile nucleon, then on the second pass one includes 
for each target nucleon so chosen, all of those nucleons from the projectile approaching it within the same transverse 
distance. This totality of mutually colliding nucleons, at more or less equal transverse displacements, constitute a 
group. The procedure partitions target and projectile nucleons into a set of disjoint interacting groups as well as a set 
of non-interacting spectators in a manner depending on the overall geometry of the A-l-B collision. Clearly the largest 
groups in P-l-A will, in this way, be formed for small impact parameters b; while for the most peripheral collisions the 
groups will almost always consist of only one colliding NN pair. Similar conclusions hold in the case of A+B collisions. 

In stage II of the cascade we treat the entities which rescatter as prehadrons. These prehadrons, both baryonic 
or mesonic in type, are not the physical hadron resonances or stable particles appearing in the particle data tables, 
which materialise after hadronisation. Importantly prehadrons are allowed to interact starting at early times, after 
a short production time j27\ . nominally the target-projectile crossing time Tab '^ Rab/i- The mesonic prehadrons 
are imagined to have {qq) quark content and their interactions are akin to the dipole interactions included in models 
relying more closely on explicit QCD |27Ll28l| . but are treated here as colourless objects. 

Some theoretical evidence for the existence of comparable colourless structures is given by Shuryak and Zahed |22| 
and by certain lattice gauge studies [33 ■ In these latter works a basis is established for the persistence of loosely 
bound or resonant hadrons above the QCD critical temperature Tc to T ~ (1.5 — 2.0) x Tj.. This implies a persistence 
to much higher transverse energy densities p{E) ^ (1.5 — 2.0)'*/5c, hence to the early stages of a RHIC collision. 
Accordingly we have incorporated into stage II hadron sized cross-sections for the interactions of these prehadrons, 
altough early on it may in fact be difficult to distinguish their colour content. Such larger cross-sections indeed appear 
to be necessary for the explanation of the apparently large elliptical flow parameter found in measurements [sj, jsj • 

The prehadrons when mesonic may consist of a spatially close, loosely correlated quark and anti-quark pair, are given 
a mass spectrum between ttt,^ and 1 GeV, with correspondingly higher upper and lower limits allowed for prehadrons 
including strange quarks. The Monte-Carlo selection of masses is then governed by a Gaussian distribution, 

P{m) = exp{—{m — mo)^/w^), (3) 



with mo a selected center for the prehadron mass distribution and w = too/4 the width. The non-strange mesonic 
prehadrons is taken at TOq ~ 500 MeV, and for strange at TOq ~ 650 MeV. SmaU changes in itiq and w have httle 
effect since the code is constrained to fit hadron-hadron data. 

Too high an upper hmit for rriQ would destroy the soft nature expected for most prehadron interactions when they 
finally decay into 'stable' mesons. The same proviso is in place for prebaryons which are restricted to a mass spectrum 
from niN to 2 GcV. However, in the present calculations the prebaryons are for simplicity taken just to be the normal 
baryons. The mesonic prehadrons have isospin structure corresponding to p, uj, or K* , while the baryons range across 
the octet and decuplet. 

Creating these intermediate degrees of freedom at the end of stage I simply allows the original nucleons to distribute 
their initial energy-momentum across a larger basis of states or Fock space, just as is done in string models, or for 
that matter in partonic cascade models. Eventually, of course, these intermediate objects decay into physical hadrons 
and for that purpose we assign a uniform decay width ^ Tf, which then plays the role of a hadronisation or formation 
time. 

B. Groups 

Energy loss and multiplicity in each group of nucleons is estimated from the straight line collision history. To 
repeat, transverse momentum of prebaryons is assigned by a random walk having a number of steps equal to the 
number of collisions suffered. The multiplicity of mesonic prehadrons cannot be similarly directly estimated from the 
number of NN collisions in a group. We argue [33 that only spatial densities of generic prehadrons f5l,i6| below some 
maximum are allowable, viz. the prehadrons must not overlap spatially at the beginning of stage II of the cascade. 
The KNO scaled multiplicity distributions, present in our NN modeling are sufficiently long-tailed that imposing 
such a restriction on overall multiplicity can for larger nuclei affect results even in P-|-A or D-l-A systems. In earlier 
work [^ Q the centrality dependence of dN/drj distributions for RHIC energy Au-|-Au collisions was well described 
with such a density limitation on the prehadrons. 

Importantly, the cross-sections in prehadronic collisions were assumed to be the same size as hadronic, e. g. meson- 
baryon or meson-meson etc. , at the same center of mass energy, thus introducing no additional free parameters 
into the model. Where the latter cross-sections or their energy dependences are inadequately known we employed 
straightforward quark counting to estimate the scale. In both SPS Pb-fPb and RHIC Au-|-Au events at several 
energies it was sufficient to impose this constraint at a single energy. The inherent energy dependence in the KNO- 
scaled multiplicities of the NN inputs and the geometry then take over. 

C. High Transverse Momenta 

One question which has yet to be addressed concerns the high p± tails included in our calculations. In principle, 
LUCIFER is applicable to soft processes i. e. at low transverse momentum. Where the cutoff in p± occurs is not 
readily apparent. In any case we can include high p± meson events through inclusion in the basic hadron-hadron 
interaction which is of course an input rather than a result of our simulation. Thus in Fig(2) we display the NSD 
{l/2TTpj_){(PN'''^°'^3'"^/dpj_ dif) from UAl (2j|. One can use a single exponential together with a power-law tail in p_L, 
or alternatively two exponentials, to achieve a fit of the output in PP to UAl y^=200 GeV data. A sampling function 
of the form 

/ = p^{aexp{-p^/w) + 6/((l + (r/a)^)), (4) 

gives a satisfactory fit to the PP data in the Monte-Carlo. 

This PP pj^ spectrum, inserted into the code, is then applied to the meson p±^ distribution in D-|-Au. No correction 
is made for possible energy loss in stage I, an assumption parallel to that made by the BRAHMS and all other RHIC 
experiments, in analysing p±^ spectra and multiplicities irrespective of low or high values. Since we impose energy- 
momentum conservation in each group, a high p±^ particle having say, several GeV/c of transverse momentum, must 
be accompanied in the opposite transverse direction by one or several compensating mesons. Such high-p^ particles 
are not exactly jets, to the extent that they did not originate in our simulation from hard parton-parton collisions, 
but they yield the same observable experimental behaviour. 



D. Initial Conditions for II 

The final operation in stage I is to set the initial conditions for the hadronic cascade in stage II. The energy- 
momentum taken from the initial baryons and shared among the produced prehadrons is established and an upper 
limit placed on the production multiplicity of prehadrons and normal hadrons. A final accounting of energy sharing 
is carried out through an overall 4- momentum conservation requirement. We emphasize that this is carried out 
separately within each group of interacting nucleons. 

The spatial positioning of the particles at this time could be accomplished in a variety of ways. We have chosen to 
place the prehadrons in each group inside a cylinder, initially having the longitudinal size of the nucleus, for a P+A 
collision, and having the longitudinal size of the interaction region at time Tab in an A+B collision, then allowing 
the cylinder to evolve freely according to the longitudinal momentum distributions, for a fixed time r/, defined in 
the rest frame of each group. At the end of r/ the multiplicity of the prehadrons is limited so that, if given normal 
hadronic sizes ~ (47r/3)(0.8)'^ fm'^, they do not overlap within the cylinder. 

Up to this point longitudinal boost invariance is completely preserved, since stage I is carried out using straight line 
paths. The technique of defining the evolution time in the group rest frame is essential to minimizing residual frame 
dependence which inevitably arises in any cascade, hadronic or partonic,when transverse momentum is considered 
due to the finite size of the colliding objects implied by their non-zero interaction cross-sections. 

III. STAGE II: FINAL STATE CASCADE 

Stage II is as stated a straightforward cascade in which the prehadronic resonances interact and decay as do any 
normal hadrons present or produced during this cascade. Appreciable energy having being finally transferred to the 
produced particles these 'final state' interactions occur at considerably lower energy than the initial nucleon-nucleon 
collisions of stage I. As pointed out, during stage II the interaction and decay of both prehadrons and hadrons 
is allowed. In the case of D+Au, although less abundant than with a more massive projectile, these final state 
interactions are as we will see, nevertheless still of some relevance. 

We are then in a position to present results for D+Au collisions. These appear in Fig(l), as previously referred 
to, and Figs(2~6), some of which are comparisons with the measurements of both BRAHMS and PHOBOS |2j3- -'-'^ 
fact the plot of experimental data in Fi gfl) is from PHOBOS [J- This PHOBOS reference also exhibits comparisons 
with several theoretical calculations |9l l35L 1371 l3a | . Two of these references 37, 38], describe much of the measured 
77-distribution at negative 77 near the target, while one |3q apparently accounts for the extreme backward tail; this is 
a subject to which we will return. 

The initial conditions created to start the final cascade could have perhaps been arrived at through more traditional, 
perhaps partonic, means. The second stage would then still proceed as it does here. We reiterate that our purpose 
has been to understand to what extent the results seen in Figures (1-7) are affected by stages I and II separately. 
i. e. do they arise from initial or from final state interactions. 

IV. RESULTS: COMPARISON WITH DATA 

Fig(3) contains the simulated charged transverse momenta spectra for D+Au at rj — (0, 3.2) alongside the UAl 
data. The many orders of magnitude fall in p± densities with rapidity is apparent. Aside from low p± the D-|-Au 
curves for increasing p± appear roughly parallel to PP; small but interesting deviations show up when the ratios 
previously defined are displayed. Additionally, direct integration of the spectra indicates that > 90% of the charged 
rapidity densities result from p± < 0.7 GeV/c. Having built in no energy loss effects on these p± distributions in 
the initial state, that a similar fall off obtains in both D-|-Au and PP is all but preordained, and the overall ratio 
between rj = 0.0 and rj — 3.2 seems to be driven completely by low p±^ physics. In fact the multiplicity choice at given 
pseudorapidity and transverse momentum is only mildly influenced by p± dependence aside from that already present 
in the PP input. 

In Fig(4) the calculated LUCIFER ratios R[DAu/PP] are plotted alongside those for BRAHMS Q at both 77 = 
and 77 = 3.2. The theoretical results are obtained using NcoU — 7.0 rather than the value closer to 7.2 employed 
by BRAHMS. Our calculation of the average number of collisions in minimum-bias D-t-Au, defined as 6 < 16 fm, is 
approximately 7.0. The Cronin affect is evident in the calculated rj ~ spectrum, less so for 77 = 3.2. This is not 
unexpected. 



A. Jet Suppression 

A very interesting result is obtained by turning off the final cascade, i. e. stage II. Then the prehadrons produced 
in stage I evolve or decay into stable particles after the time tj ^^ l/(r) and do not otherwise suffer interaction. 
This situation is described in Fig(5), where it is clear that the magnitude of the Cronin enhancement of dN/dp± 
is considerably magnified. This enhancement is then very much a creature of I, i. e. a product of the transverse 
momentum gained in collisions with nucleons in the target. Incidentally, Fig(5) also indicates that a compensating 
increase in transverse momentum density occurs at the lowest p±. 

One might well turn this around and declare that the final state scattering of a given prehadron with comovers 
has cut down the Cronin effect. This is an appreciable reduction which suggests the applicability of the term 'jet 
suppression', a reduction which indeed constitutes final state suppression. The change in p± spectra is considerably 
less for 77 ~ 3 where Fig(l) indicates considerately less comovers are present, and indeed the Cronin enhancement is 
less evident at the more forward ry. 

The spillover of such stage II comovers decreases with increasing distance from the target pseudorapidities. The 
Cronin rise is itself less evident at forward rj, the ratio R[DAu/PP] is flatter as a function of p±. Indeed as rj 
increases the p± spectra approach to the PP spectra. This is again easily understandable: the most peripheral 
collisions involving the least number of participants will contribute more strongly to more forward rapidities. Both 
the unrenormalised theoretical and measured D+Au curves appear to merge with PP at the largest pseudorapidities 
shown in Fig(l). One expects to see a corresponding behaviour with decreasing centrality and decreasing participant 
nucleon number. 

B. Centrality Dependence of R[DAu/PP] 

A similar theme then is repeated in Fig(6) where the p± spectra for two quite different degrees of centrality differ 
markedly; the disparate centralities are defined by 6 < 4 fm and b = 8 fm. The 77 dependence exhibited for 6 = 8 
fm, a clearly peripheral geometry, is very muted with both rj = and rj = S showing a strong resemblance to the 
PP ?7 = data. The more central choice, 6 < 4, is subject to quite strong pseudo-rapidity variation. This explains 
at least qualitatively, the behaviour of the second ratio BRAHMS focuses on, i. e. Rep. Q]- The crossover of this 
ratio with centrality as a function of ry is explained by the much flatter dependence in rj discussed here. The differing 
number Ncoii then must be invoked to complete the picture but it plays only a passive role, the flattening of the 
p^-distribution already evident in Fig(6) for widely differing centralities or impact parameter plays the essential role. 

V. CONCLUSIONS 

It is hard to conclude definitively from what is presented here that the gluon saturation |3J, |33 and attendant 
colour-glass-condensate interpretation 183, [Sg of the BRAHMS data is not a more fundamental explanation of the 
measurements discussed here. Certainly the low- a: basis for this modeling is related to the increasing rj picture 
presented here, and perhaps the gluon saturation aspect of that approach is mirrored in, and underpins, the prehadron 
multiplicity limitation employed above. 

It would seem however that the direct attempt at a PQCD explanation of this behaviour must claim that, at the 
very least, all soft mesons are produced in essentially hard collisions. The presentation here provides an interesting 
case for relying on essentially soft, low p±, processes to produce the major features of the BRAHMS data. True 
enough, the high p± tails in distributions are merely tacked on in our approach, but legitimately so by using the 
PP data as input to the nucleus-nucleus cascade. The PP dN/dp± variation with p± to a large extent drives the 
multiplicity generation in D-|-Au, altering only slightly the hard ry-dependent ratio from the soft. One only need add 
the assumption that the p± tails in D+Au do not exhibit any drastic non-monotonic behaviour. 

The nascent appearance of appreciable high p± suppression, especially in our rj—O spectrum seen in Fig. 5, suggests 
that enhanced suppression will occur in a full Au-)-Au collision. Whether the complete, or an appreciable fraction of, 
jet suppression seen in RHIC experiments can be explained by final state interactions remains to be established. We 
note that C. Greiner and coauthors !2l| and Gassing |22| have commented forcefully on precisely this point, also in a 
hadron-based cascade setting. 

We indicated we would return to the target pseudorapidity region. It is of direct relevance to do this for discussion 
of specifics in D-|-Au but also for the implications for Au-|-Au and other complex systems. For the deuteron projectile 
PHOBOS data not only extends further backward than other experiments or calculations, but also exhibits a feature, 
perhaps a shoulder, near where our calculations exhibit a peak in the charged baryons (see Fig(I)). It would clearly 
be of some import to have particle identification in the present measurements of rj and p± distributions. Considering 



the D+Au system, one notes that transverse momentum distributions near the target, or further back in rj, are 
significantly softer, again possibly anticipating high p± suppression to be associated with the symmetric massive ion 
collisions. 
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FIG. 1: D+Au Pseudorapidity spectra: Direct comparison of PHOBOS minimum bias data with LUCIFER simulation, the 
latter for blelG fm. Two calculations are shown, slightly differently normalized, one for which the average prehadron decays 
into 3.15 observed mesons and one for which this number is 3.8. Both are only slightly above that expected for pure NN 
production, indicating some calculated events exceed the multiplicity constraint discussed in the text. The absence of collision 
number divisors is instructive, revealing both the considerable production of final state mesons at ?7=0, in excess of PP, and 
the apparent asymptotic approach of both data and calculation to PP at the largest observed rj. 
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FIG. 2: Transverse Momentum Spectra; UAl vs LUCIFER. The basic input of both low and high pt for simulations obtained 
from a fit to the SPS data from the UAl Collaboration l23l. 
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FIG. 3: The simulated charged transverse momenta spectra for D + Au aX r] = (0, 3.0) alongside the UAl data. Fits to 
the LUCIFER calculations, used to interpolate the simulation are shown. These are made with the combination of a single 
exponential at low pr and a power law at higher values. 
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FIG. 4: Minimum Bias RDAu/PP for r;=(0,3.2): The BRAHMS results are compared to the collision number-normalized 
calculations. The latter are obtained using the results in Fig. 3 with Ncoii=7-0, compared to the BRAHMS choice 7.2 ± 0.3. 
The presence of a Cronin effect is clear, with however a flatter pr dependence obtaining for the larger tj. 
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FIG. 5: EfTect of the Final State Hadronic Cascade. The Cronin enhancement of d N/dr)dp±_ for 77=0 is markedly reduced 
by the second stage of LUCIFER. For larger p± one might refer to this phenomenon as final state 'jet' suppression. The 
diminuition is muted for larger 77. For more massive ion-ion collisions one can expect a considerably greater reduction. 
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FIG. 6: Centrality Dependence of RDAu/PP. The two sets of r;=(0,3.2) transverse momentum distributions for one central 
6 < 4.0 and one more peripheral impact parameter fe=8.0 are displayed against the UAl PP data. Evidently the rj dependence 
is strong for the central choice and virtually vanishing for the peripheral collision, with the latter distributions closely matching 
UAl. 



